Rotational Symmetry Breaking in Sodium Doped Cuprates 

Yan Chen\ T. M. Rice^-^, and F. C. Zhang^ 

O ' ^Department of Physics and Center of Theoretical and Computational Physics, 

O 

CN ' The University of Hong Kong, Pokfulam Road, Hong Kong, China 

U '. 

Qh' Institut fiir Theoretische Physik, ETH-Ziirich, CH-8093 Switzerland 



< 



o 



o 
o 



% 



(Dated: February 6, 2008) 

Abstract 

For reasonable parameters a hole bound to a Na+ acceptor in Ca2_a;NaxCu02Cl2 has a doubly 
degenerate ground state whose components can be represented as states with even (odd) reflection 



c/3 ■ symmetry around the x{y) -axes. The conductance pattern for one state is anisotropic as the tip 

-4— > . 

Cd . of a tunneling microscope scans above the Cu-O-Cu bonds along the x(y)-axes. This anisotropy 

is pronounced at lower voltages but is reduced at higher voltages. Qualitative agreement with 

Ch ' recent experiments leads us to propose this effect as an explanation of the broken local rotational 

^ ■ symmetry. 
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The recent atomically resolved scanning tunneling microscopy (STM) studies by Davis 
and collaborators [l|, |^ on strongly underdoped Ca2-a;Na2;Cu02Cl2 revealed a surprisingly 
complex pattern with the square symmetry of the lattice broken on a local scale. The STM 
data were analyzed to obtain the local hole density of Cu-site using the method proposed 
by Randeria, Sensarma, Trivedi and Zhang |3] and also a related method proposed by 
Anderson j^. In these methods, the differential conductance signal is integrated from the 
chemical potential to a substantial voltage cutoff. Randeria et al. [3] pointed out that the 
ratio between its positive voltage (electron injection) and negative voltage (hole injection) 
signals measured the local hole density of the Cu-site independent of the strength of the 
troublesome tunneling matrix elements. 

In Ca2_i.Naa;Cu02Cl2 the acceptor Na"''-ions substitute for Ca^+-ions and sit in the center 
of a square of four Cu sites both above and below the outermost Cu02 layer. The attractive 
potential for the doped holes generated by a Na"*" acceptor does not break the local square 
symmetry. The local square lattice symmetry broken observed in the STM data is most 
pronounced at the in-plane oxygen sites and for the lower voltage cutoff. This has raised 
an interesting question about the origin of the broken local square symmetry in the STM 
experiments. 

In this letter we address the origin of this broken local symmetry. We show that such 
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FIG. 1: Schematic crystal structure of Ca2_xNaxCu02Cl2 around a Ca^"*" ion or a doped Na+ ion. 
The sign represents the phase of the Cu- and Cl-orbital wavefunction. 



broken symmetry states appear in the case of a hole confined to a cluster of sites centered 
at a Na"*"-acceptor. Further this broken local symmetry shows up in the STM mapping 
primarily as a modulation pattern above the 0-sites in the outermost Cu02-layer. 

We start by considering the STM tunneling process in detail in Ca2_xNa^Cu02Cl2. As 
pointed out by Rice and Tsunetsugu J5[ , the STM tip will couple primarily to the Spz-states 
of the outermost CI layer but the matrix element for an electron to tunnel directly into a 
hole located on the Cu02 square directly underneath the Cl-ion vanishes by symmetry. The 
holes move in a band of singlets in a state with d^2_j^2 -symmetry centered on a Cu-site but 
strongly hybridized with the four nearest-neighboring (n.n.) O sites that surround each Cu. 
The 3p2-Cl states however hybridize with the hole states centered on the four n.n. Cu sites. 
Each of these states hybridizes either through the direct overlap of the Sp^-Cl with 2p^.(j^)-0 
orbitals and also through the 4s-Cu orbital on the site directly below each CI. Let I^E'J'^) 
denote the single hole ground state with energy E^^, we can write the amplitude to inject 
an electron in a Sp^-Cl hole wavefunction as a superposition in hole-states centered on the 
four n.n. Cu-site. 
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where p ^ (cl ^ ) are creation operators for electrons in the 3p^-Cl orbital at site i (planar 

coordinate) and the d-p hybridized orbital centered on the Cu-site at i + r in the outermost 

Cu02 layer respectively, r = ±x, ±y is a planar vector connecting n.n. Cu sites, (i, Cl\i + 

r, Cu) denotes the overlap between a Sp^-Cl and the d-p hybridized orbital centered at a 

n.n. Cu site. Assuming the square lattice is four-fold rotational invariant, the overlaps are 

independent of r except for their sign, which leads to the final result in Eqn. (1). For d^^-y^ 

symmetry of the Cu-orbital, we have (—1)*^ = —1 for r = ±x, and (—1)^^ = +1 for r = ±y. 

Note that only the relative phase is important. 

Following the theory of STM tunneling processes developed by Tersoff and Hamann J6[, 

we can write the differential conductance at voltage V at r, the center of curvature of the 

tip, 

dl{f) 



dV 



ocEK^l4,JOlV-i.o-.- (2) 



where al^ is the electron creation operator at position f, \m) are eigenstates of the half-filled 
system with energy Em, and uj = eV. When the tip is scanned from above the CI site at i 
to a neighboring site, i + r', (planar coordinates), we have 

The integrated current at r up to a positive voltage V is then, 

(4) 
where A is a constant, and G is a step function. The variation of the tunneling current as 
one scans between two n.n. CI sites will therefore be sensitive to the relative phase to inject 
electrons on n.n. Cu sites in {"^l^)- To this end we examine the lowest energy single-hole 
states for clusters up to 16 sites including an attractive potential for the center square. 

The two-dimensional single-band t-t'-J model with the on-site impurity potential on four 
Cu sites of the center square around the Na"*" dopant is defined by the Hamiltonian, 

i,<y {i,j)a 

-t' J2 {clc,^ + h.c.) + jJ2Si-Sj. (5) 

A constraint of no double electron occupation is implied. {i,j) and {{i,j)) refer to n.n. and 
next n.n. sites i and j. The onsite potential of the four n.n. Cu sites around the sodium 
dopant is e^ = e^ while e^ = for the rest of sites. The single hole impurity state of the 
Hamiltonian (4) at t' = was studied previously by using exact diagonalization method for 
small clusters by Von Szczepanski et al. '7|, Rabe and Bhatt \^, and Gooding [9|. The 
ground state was found to be orbitally doubly degenerate (5* = 1/2) in certain parameter 
regime [8|, |9|, |lO| . When the next n.n. hopping integral t' is included, our results below show 
that the doubly degenerate S=l/2 states is lowest in energy within a reasonable parameter 
range. The symmetry of the lowest energy state however is sensitive to boundary condition 
and parameters. In our calculations, we choose t = 1 and J = 0.3, suitable for the cuprates, 
and study the ground state as a function of t' and e^. We consider the clusters of 12- and 
16-sites as illustrated in Fig. 2. 

Exact diagonalization calculations of the t-t'-J model on small clusters have been re- 
ported lul, 112, llSl llJ, ll5|, llg . For a 16- (4x4) site cluster with periodic boundary condition 
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FIG. 2: The configuration of 12-site (a) and 16-site cluster (b). The dashed line corresponds to 
the reflection axis for Px(y)- The red cross marks the tip position on the middle point of two n.n. 
Cu sites along the x- and y-axes. In (b), eight independent bonds are labeled with numbers. 

(PBC), the ground state of a single hole has a four-fold symmetry with hole momentum 
(±7r/2, ±vr/2) if t' < 0, and a two-fold symmetry with hole momentum (vr, 0) or (0, vr) if 
t'>0. 

The presence of an impurity potential on the lightly doped t-J model affects not only 
the local charge and spin distributions but also the symmetry of the ground state wave- 
function jll, y, |9| . In this study, we focus on the reflection symmetries of a two-dimensional 
square lattice with respect to x- and y-axes passing through the sodium dopant {P^ and Py 
respectively) and on the parity PxPy Since [Px{y),H] = 0, we may classify states according 
to the quantum numbers of Px, Py We denote the state with [Px = +l,Py = +1) as 
phase (I), and doubly degenerate state (+1,-1), and (—1,+!) as phase (II), and (—1,-1) 
as phase (III). In Fig. 3(a) and 3(b), we show the ground state phase diagram in terms 
of its reflection symmetry obtained by exact diagonalization. As we can see, there is a 
regime with the doubly degenerate ground state. Previous band structure calculation [iTj 
and angle resolved photoemission experiments |l8| both indicate a rather weak t' ~ —0.1 for 
Ca2-a.Naa;Cu02Cl2. Thus we expect the existence of doubly degenerate ground state in this 
doped material. We note that the four low-lying states with distinct reflection symmetries 
are very close in energy. Small perturbation of lattice geometry, boundary condition, t' or e^ 
may change symmetry of the ground state. In Fig. 3(c), we show the energy level evolution 
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FIG. 3: Ground state phase diagram in the parameter space of the next n.n. hopping integral i! 
and the sodium impurity potential e^ in (a) 12-site cluster (black solid line) and 16-site cluster (red 
dashed line) with open boundary condition, (b) 16-site cluster with PBC. Regime D: the doubly 
degenerate state; regime S: the singly degenerate state. The thick black line in (b) corresponds 
to the four- fold symmetric states at e^ = 0. (c) Energy levels as functions of e^ at t' = —0.12 as 
indicated by the dashed green line in (b). Ei denotes the energy in phase I. The dashed line is for 
doubly degenerate states while the dotted line is for phase III. 

as e^ changes for a fixed value of t' = —0.12. 

We next examine the effect of the rotational symmetry broken ground state to the STM 
measurements. We consider the system to be in one of the two degenerate states, say in the 
state of Px = l,Py = —1. The broken symmetry may be caused by quadrupole interaction 
of two single-hole states or by other couplings. We remark that the local density of states at 
the Cu-sites remains rotational symmetric in this state. However, because of the interference 
of different Cu-sites to the integrated differential conductance in Eqn. (4), the STM tip at 
a position between two Cu-sites will break the rotational symmetry. To be specific, we 
consider the tip position above the two geometrically symmetric points at halfway between 
the Cu-site i and its neighboring sites i + ct or i + y. The sodium impurity is located at 
i + x/2 + y/2. Note that {f\i, CI) = (r|i + r', CI) for the mid point r, we obtain the expression 
for the integrated differential conductance P^y\iu) = I{^{^),uj) at position i + 



f(|) with 



the cut-off energy a;, 

a 

if'H-^) = Ao 5^ Km| $:(-l)*^^l4^.^ + ii.(,)+.>)l^o')r0(^ -Em + El\ (6) 

m f 

where la is the spin-dependent integrated differential conductance, and the Cu-site index i 
has been dropped for simphcity. The conductance can be calculated in a small cluster. We 
diagonalize the hamiltonian H in Eqn. (5). exactly and find the ground state of the single 
hole and also all the eigenstates \m) and the corresponding energies E^ at the half- filling, 
which allows us to calculate all the matrix elements in Eqn. (6). Note that I^{uj) = P{uj) 
by symmetry if |\E'J'^) is non-degenerate. In the parameter regime where I^E'o'^) is doubly 
degenerate (see Fig. 2), I^ and /^ can be different. In Fig. 4(a), we plot I^{u) and P{u) 
as functions of the cutoff energy u obtained from the 12-site cluster calculations, where the 
hole groundstate has symmetry P^ = +1, Py = —1 and has a spin 82=1/2. The asymmetry 
of I between the two tip scans is apparent and is more pronounced at a lower energy cutoff 
and becomes weaker at a high energy cutoff. To see the asymmetry more clearly, we also 
plot the ratio of the conductances at the two tip points in the same figure, 

v{uj) = P{u)/P{u;). (7) 

It is interesting to note the strong cutoff energy dependence of the asymmetry in the 
conductance. For example, rj ^ 0.3 for uj ^ J, and rj ^ 0.85 for uj = AJ. At uj —>■ 00, 
1] ^ 1.04, which shows a very weak asymmetry. Our theoretical results are consistent 
with the recent STM observation J], where the strong asymmetry has been found at a 
lower energy 150 meV (~ J), and the asymmetry becomes much weaker or unobservable 
at a high energy 600 meV (~ AJ). In Fig. 4(b) and 4(c), we show the spin- dependent 
conductances /o-(r'/2,Lo') as functions of uj in different spin channels. Our theory predicts 
that the asymmetry is spin-dependent. At high energy cutoff limit, we find the ratio P/I^ ~ 
1.23 in the spin-up channel and P/If ~ 0.88 in the spin-down channel. The asymmetry 
compensates each other in the two different spin channels, which gives a weaker asymmetry 
in the total ratio P/P ~ 1.04. We remark that rapid spin flip processes around the sodium 
impurity may make the spin-resolved STM experiments more difficult to observe. However, 
at very underdoped samples, the spin may be frozen around the sodium impurity, which 
should provide the possibility for the spin-resolved STM experiments. 
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FIG. 4: The total integrated differential conductance on two cross (X) sites in Fig. 2a along the x- 
(solid line) and along the y-axis (dashed line) as functions of the cut-off energy lv together with the 
anisotropy ratio rj (thick solid line). The results for spin- up channel (b), and spin-down channel 
(c). In the calculations (b), and (c), the single hole state has a spin 82=1/2. 



To better understand our results, we analyze the groundstate and the conductance in a 
four-site cluster system. The 4-site cluster of the Hubbard model with a single hole was 
previously studied by Altman and Auerbach 19|. They have found that for not very large 
on-site Coulomb repulsion U, the ground state is two-fold degenerate with spin-1/2. We 
here study the single hole ground state of the Hamiltonian Eqn. (5), and focus on the 
spin-1/2 sector, which is relevant for not too small ratio of J/t. Note that the impurity 
potential does not play any role in the 4-site cluster. The ground state is doubly degenerate 
with the energy given by Eg = -J/2 - ^/WTWTJJW^ with tan^i = Vst/{J + t' + 



Eg). The two degenerate states can be classified by symmetry P^ = +1,-Pj 



-1 and 



Px = —l,Py = +1. Assuming the former to be ground state, we find that apart from 
an overall prefactor, lUuj) = |cos^70(ct; — J), If{uj) = |cos^70(ct; — J), I^{uj) = and 
Jf (cu) = cos^ 7G(c(j)+sin^ 7G(c(j— 2J). Due to the nature of the many body wavefunction, the 
conductance in Eqn. (6) may involve either destructive or constructive interference among 
different components. This effect is most drastic in H with an exact cancelation between 



them. Summing over spins, we have, P{u!) 



JX _|_ JX 



2cos^-fe{uo - J) while iy{uo) 



bond index 
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0.342 


0.328 


0.122 


0.129 


0.020 


0.014 


0.012 


0.016 


-{S.S,) 


0.167 


0.074 


0.161 


0.167 


0.329 


0.336 


0.347 


0.334 



TABLE I: The hopping integral and spin-spin correlation function for various bonds in a 16-site 
cluster with PBC for t' = —0.1 and e^ = 1.0. The bond indices are labeled in Fig. 2(b). 

cos^ ■yQ{uj) + sin^70(c<j — 2J). It is obvious that rj^uj) becomes infinite for u G (0, J), ri{uj) = 
1/2 for uj G (J, 2 J) while //(tu) = 2 cos^ 7 for cj > 2 J. At the limit of the large bias voltage, we 
have 1] = 2 cos^ 7. For t'/t = —0.1, rj ~ 0.97. Therefore, the asymmetry becomes rather weak 
with the increase of bias voltage. The ratio of the integrated conductances in the ground 
state with symmetry P^ = —l,Py = +1 can be obtained by the interchange symmetry of x 
and y, and is given by l/rj. We have also calculated the integrated conductance for the hole 
injection tunneling, and found also rotational asymmetric. The ratio of the conductances 
between the electron injected and hole injected is also asymmetric. 

In the rotational symmetry broken state we considered above, the local charge density at 
each of four n.n. Cu sites around sodium dopant is uniform, but the n.n. Cu-Cu bondings 
are not equal. For instance, in the ground state with Pj. = +1, Py = — 1 of the 16-site cluster, 
we show the expectation value of hopping integral (CJCj) as well as spin-spin correlation 
function (SiSj) for various bonds, depicted in Table I. The hopping integral along bond 1 is 
slightly stronger than that of bond 2 while the spin-spin correlation function exhibits a much 
stronger asymmetry. The impurity may lead to localization of holes around this impurity 
which may result in a larger hopping integral and weaker spin-spin correlation function 
around this impurity. Close to the lattice boundary, since the hole density becomes very 
small for large impurity scattering strength, the spin-spin correlation function may approach 
to -0.34 and hopping integral vanishes. Our results are consistent with previous calculations. 
The strong variations of hopping integral and spin-spin correlation around the impurity may 
lead to the formation of local lattice distortion and the appearance of spin-lattice polaron 
due to the electron-phonon coupling [20|, which may induce further rotational asymmetry. 
We speculate that this could be one of the reasons for the inhomogeneity in underdoped 
Ca2„,.Na,Cu02Cl2. 
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